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Turbulence induced by the ion temperature gradient (ITG) is investigated in the helical and axisymmetric 
plasma states of a reversed field pinch device by means of gyrokinetic calculations. The two magnetic con¬ 
figurations are systematically compared, both linearly and nonlinearly, in order to evaluate the impact of the 
geometry on the instability and its ensuing transport, as well as on the production of zonal flows. Despite its 
enhanced confinement, the high-current helical state demonstrates a lower ITG stability threshold compared 
to the axisymmetric state, and ITG turbulence is expected to become an important contributor to the total 
heat transport. 


I. INTRODUCTION 

A long standing question in the context of magnetic 
fusion is the impact of 3D shaping on the fundamental 
properties of the plasma. The two most studied toroidal 
configurations, stellarators and tokamaks, manifest in 
general different magnetohydrodynamic (MHD) stability 
and neoclassical confinement. Comparisons of this sort 
have already been addressed^, and the underlying reasons 
seem to be well understood. Conversely, the difference in 
the behavior of turbulent transport is a topic which re¬ 
mains to date largely unanswered: while the tokamak 
line has been extensively investigated in this regard dur¬ 
ing the last decades, it is only recently that the stellara- 
tor community has started taking sophisticated steps in 
this direction, due to the difficulty in dealing with helical 
fields^. 

In this paper we exploit the unique flexibility offered by 
the reversed held pinch (RFP) device, thanks to its abil¬ 
ity to produce both helical and axisymmetric plasmas in 
the course of a single discharge. The RFP is nominally an 
axisymmetric configuration, which at low plasma current 
is characterized by a wide spectrum of resonant MHD 
modes maintaining an overall (quasi)axisymmetry of the 
plasma. Interestingly though, the RFP plasma at high 
currents experiences a transition to a helical state, named 
single helicity^. Depending on the current intensity, the 
duration of the helical state can be long enough to reach 
an equilibrium state, which can be described numerically 
by special codes like VMEC—. 

Despite the fact that the stability of the ion temper¬ 
ature gradient (ITG) nrode&£ has already been studied 
in the RFP, so far only axisymmetric equilibrium models 
have been taken into account—The main conclusion 
from these geometrically simplified investigations is that 
the ITG stability threshold in the RFP is larger than 
in the tokamak, typically by a factor R/a (with R the 
major radius and a the minor radius of the torus). An 
explanation comes from a detailed analysis of the paral¬ 
lel dynamics, showing a relevant Landau damping of the 
mode due to the short field connection length— Nonlin¬ 
ear simulations of ITG turbulence have been performed 


with and without impurities^, showing a relatively low 
ion heat transport and a significant Dimits shift. 

In this work we aim at revisiting these findings, consid¬ 
ering what the introduction of a helical deformation may 
cause to ITG turbulence. To do this, we rely on the 
massively-parallcl Eulerian gyrokinetic code GENE—, 
applied to the VMEC helical/axisymmetric RFP equi¬ 
libria with the aid of the code GIST—. After introducing 
the MHD equilibria of the RFP in Sec. El we present 
the comparison of the respective ITG modes in Sec. Ell 
to proceed with the ITG turbulence and the behavior 
of zonal flows in axisymmetric and helical systems in 
Sec. ED We conclude with a short discussion in Sec. El 
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FIG. 1. Safety factor profiles (top) and cross-sections at 
toroidal angle (f> = 0 (bottom) for the helical (left) and ax¬ 
isymmetric (right) VMEC equilibrium reconstructions. The 
grey-shaded areas depict the selected radial domain for the 
gyrokinetic investigation. 
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II. AXISYMMETRIC AND HELICAL EQUILIBRIA 

The RFP configuration is characterized by a low safety 
factor profile, q < 0.1 in the core, which further decreases 
in the outer region, reaching slightly negative values at 
the very edge. In case the plasma can be assumed ax- 
isymmetric, the connection length of the field turns out 
to be L c = 2nR(q 2 + e 2 ) 1 / 2 , with e = r/R, which is 
much lower than in a tokamak; in particular L c ~ 2™ 
for r ~ a, where the field is purely poloidal. 

In RFX-mod, pushing large plasma currents (> 1 MA) 
makes the plasma undergo a transition to a helical state, 
named single helieity, with a single tearing mode satu¬ 
rating well above the others. Axisymmetric and helical 
RFP plasmas have rather different features in terms of 
overall transport properties, MHD dynamics, impurity 
behaviour—. We mention here only the sharp reduction 
of magnetic field line stochasticity in the helical states, 
with the occurrence of transport barriers, at least in the 
electron heat channel. For this reason, we expect tur¬ 
bulence (at least temperature-gradient driven) to play a 
potential role in such states. 

From the point of view of MHD, the most relevant pa¬ 
rameter differentiating the axisymmetric from the single- 
helicity state is the q profile. Axisymmetric states have 
a monotonically decreasing q as a function of the radius; 
for the helical equilibria the q profile has a maximum, 
which statistically corresponds to the maximum electron 
temperature gradient—^. There has been a considerable 
effort of the RFX-mod team to describe such equilibria 
by means of the equilibrium code VMEC, and more re¬ 
cently by means of V3FI T 17 i 18 ; which uses VMEC as a 
solver to determine the equilibrium that best matches the 
experimental data. In particular, VMEC has been modi¬ 
fied to work in the RFP: in a reversed field configuration 
the toroidal flux ip t cannot be used as a radial coordinate 
in this configuration, as it is not monotonic. For this rea¬ 
son a new version of the code has been released, which 
uses the poloidal flux ip as the radial coordinate. 

However, since the magnetic field representation in 
GIST makes use of the toroidal flux as a radial coordi¬ 


nate, we prefer to slightly modify the q profile in the very 
edge, keeping it slightly positive. This approximation 
does not influence the q profile and other equilibrium- 
related quantities in the radial region where we are going 
to perform our study, i.e., the mid-radius region where 
transport barriers emerge. Thus, in the following, the 
toroidal flux normalized at the last closed magnetic sur¬ 
face, s = ipt/tpt,ids, s G [0,1], will play the role of the 
radial label. 

In Fig. [Owe show two equilibrium reconstructions. On 
the left, the VMEC reconstruction is made using V3FIT, 
therefore a minimization procedure is performed on the 
experimental data. On the right, the reconstruction is 
based on an axisymmetric q profiled, i.e., it does not 
make use of the (to, n) = (1, 7) tearing mode which drives 
the helical deformation. We prefer to use two different 
reconstructions at the same time instead of considering 
different time instants of a discharge, so as to keep most 
of the plasma parameters equal. 

The radial region where we are performing gyrokinetic 
simulations is s £ [0.3, 0.8]. The majority of them are on 
the s = 0.5 surface, where the (helical) q profile peaks 
and the magnetic shear is almost vanishing. The GIST 
code prepares a flux tube domain with the full geomet¬ 
ric description in terms of metric coefficients, curvature 




FIG. 2. A magnetic flux surface of the helical equilibrium, 
with the a = 0 (blue) and the a = 7t/7 (red) flux-tubes, 
where a = qd — <p is the field-label coordinate. The dot on 
each tube marks the origin in the parallel direction, z = 0. 


FIG. 3. Various geometric coefficients (see text for definitions) 
versus the parallel coordinate z on the surface s = 0.5, for 
the axisymmetric case and for the two stellarator-symmetric 
tubes, a = 0 and a = 7r/7, for the helical case. 
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operators, Jacobian, magnetic field strength, and parallel 
gradient- which the GENE simulation is based on. 

In principle, it is possible to obtain either a Boozer 
or a PEST coordinate representation of the flux tube. 
We will use the latter in the following, as that repre¬ 
sentation is more directly linked to the VMEC (s, i?, <j>) 
coordinates, d and <f> being the VMEC poloidal angle and 
the cylindrical toroidal angle, respectively. In particular, 
GIST makes use of the straight-field-line poloidal angle 
6 = d + A(s, i9, (f>), with A the stream-function of VMEC. 
Each point in the flux-tube domain at s = sq is described 
by a triplet ( x,y,z ), with x = s 1 / 2 (radial coordinate), 
V = ( s 1/2 /q ) (qO - 4 >)| S=S0 = (s 1 / 2 /g)| s=5o a (binormal), 
z = 0 (parallel). 

Due to the peculiar helicity of the helical states, 
(m, n) = (1,7), the overall magnetic geometry has a 
seven-fold symmetry. Two stellarator-symmetrio=* flux 
tubes are used for gyrokinetic studies, one centered in 
the outboard midplane of the (f> = 0 section, the other 
one centered in the outboard midplane at <fi = —7t/T, 
see Fig. [2] The two tubes have a = 0 and a = n/7 re¬ 
spectively. In Fig.[3]we show some relevant (normalized) 
parameters as a function of the parallel flux-tube coordi¬ 
nate the g xx = \ Vx\ 2 and g vv = | Vz /| 2 components of 
the metric tensor, the normalized magnetic field B and its 
curvature = ( By/g^K noinl +g xv K geo )/g xx (with K n0 rm 
and Kgeo the normal and geodesic components of the field 
line curvature, respectively; negative u>k corresponds to 
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FIG. 4. ITG mode growth rate (a)-(b) and real frequency (c)- 
(d) for the helical a = 0 tube (left frames) and axisymmetric 
(right frames) geometry, on the respective s = 0.5 surface. 
Here, a = ( Ti/nn ) 1 ^ 2 is the ion sound speed and a is the 
minor radius. 


unfavorable curvature causing ITG instability), the Ja¬ 
cobian J = (V* x Vy • Vz) _1 , and the local magnetic 
shear si oc = d(g xy / g xx )/dz, with the dashed lines repre¬ 
senting the values of the (global) magnetic shear, given 
by s = (sioc)z = 2 sq'(s)/q. The even-parity with respect 
to z = 0 reflects in the stcllarator symmetry of the tubes. 
The other case shown in the figure is the axisymmetric 
configuration, where, of course, a does not play any role. 


III. ITG INSTABILITIES 

In RFX-mod, a local estimate of the ion tempera¬ 
ture gradient is not available in the plasma core. Ion 
temperature profiles have been obtained recently from 
neutral-particle-analyzer data, these profiles being reli¬ 
able in the edge only, because of the increased neutral 
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FIG. 5. ITG mode growth rate 7 as a function of the bal¬ 
looning angle do for the a = 0 and a = n/7 helical flux tubes 
on the s = 0.5 surface, for binormal wavenumber k y pi = 1.2 
and ion-temperature gradient 1 /Lr t =5 (a). The maxima of 
the g xx component of the metric tensor (b) correspond to the 
minima of |ftx| (with do such that 7 (^ 0 ) is maximum, see text 
for definitions) (c), with the structure of the corresponding 
normalized eigenfunctions in (d). 
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absorption probability experienced towards the plasma 
center—. The Doppler broadening of spectral lines from 
impurities provides, too, good local estimates of T/ in 
the outer plasma only, where a ratio Tj/T e ~ 2/3 is com¬ 
monly evaluated— Therefore, we consider the electron 
temperature T e (s) as a reference profile in the core. A de¬ 
tailed study of the dynamics of the electron temperature 
barriers arising during the helical states can be found in 
Ref. [ 23 ]. In particular, the T e barriers turn out to have 
logarithmic gradients 1 /L Tb = — T~ x dT e jds% 2 ~ 2 — 3. 
Our studies on the ion temperature gradient must cover 
a larger interval. 

In the previous section we have seen how the flux-tube 
domain is properly defined. The code GENE can now 
solve the system of gyrokinetic equations to investigate 
the electrostatic ITG mode. Since we are emphasizing 
the geometrical aspects of the instability, we consider for 
simplicity adiabatic electrons. In addition, the density 
gradient 1 /L n = — n~ x dn/ds % 2 is assumed to vanish ev¬ 
erywhere. This is a good approximation in the region 
s < 0.8, as the experimental density profiles turn out to 
be essentially flat both in the axisymmetric and in the 
helical state. High density gradients do exist in the edge, 
but that domain is outside the scope of this work. Other 
assumptions are zero collisionality and plasma /?. While 
the ion collisionality v does not introduce relevant dif¬ 
ferences in the linear results (in RFX-mod experimental 
conditions, i//(u t hy/a) ~ 10“ 3 — 10 -2 , with T) ~ 2/3 T e ), 
assuming vanishing (3 can be more questionable (in RFX- 
mod, /3 ~ 1 — 2%). As is known, a finite /? stabilizes the 
ITG mode both in the tokamak and the stellarator, see, 
e.g., Ref. [24| for a recent study on this topic. In axisym¬ 
metric RFPs, the /3 suppression turns out to be moderate 
with respect to the tokamal—, with ITG modes being 
unstable until /3 ~ 10%. Whether this slow quenching is 
a feature also of a helical RFP is an open issue, left for 
a future work. 

Let us consider two cases, the helical equilibrium with 
flux tube a = 0 and the axisymmetric one, both on the 
surface s = 0.5. The other tube for the helical case 
( a = n/7) will be considered in a while. The parallel 
domain is set to cover 2 £ [—87t, +87t]. In the velocity 



FIG. 6. ITG mode growth rate as a function of the ballooning 
angle for several adjacent flux tubes. The largest growth rates 
are localized on the external ridge of the helical structure. 


space, the box size in the u|| direction is 3 Vth,i with 32 
grid points, while in the g direction 9Ti/B 0 with 8 grid 
points. These conditions are chosen both for the helical 
and the axisymmetric equilibria, even though shorter do¬ 
mains in the parallel direction are typically required for 
the axisymmetric equilibria. 

Increasing the ion temperature gradient 1 /L^ = 

— T~ l dTi/ds 1 / 2 , all the other parameters fixed, the lin¬ 
ear spectra in Fig. []] are obtained; the introduction of 
a helical deformation strongly destabilizes ITG modes. 
To understand this feature, let us recall that in the 
ballooning representation^ perturbed quantities vary as 
S = 6(6) exp (iS), where <5 is slowly varying in the parallel 
direction and S = S(ilj,a) describes the rapid variation 
across the field line. The perpendicular wave vector is 
VS 1 = k± = fc^VV’ + fcaVa, the components = d^S 
and k a = d a S being constants. Given S = k a (a — qOo ), 
with 9q the ballooning angle, it follows that fcy, = —k a q'9o 
(or, in normalized units, k x = —k y s9o, with s magnetic 
shear). Thus k\ = kl(g aa - 2 q'9 0 g^ a + q' 2 9lg ^). For 
9 0 = 0, which is the case of Fig. Q] (this choice will be jus¬ 
tified later), the wave vector is reduced to k 2 ± = k^g 0101 . 
As can be seen in Fig. [3J the metric component g aa is 
much smaller in the helical case than in the axisymmet¬ 
ric case. This implies a reduced finite Larmor radius 
(FLR) suppression of the ITG mode via the Bessel func¬ 
tion Jo(k±v±/£l) in the linear gyrokinetic equation. Even 
if the helical curvature is not so unfavorable as the ax¬ 
isymmetric one -about one half where the mode peaks- 
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FIG. 7. ITG growth rate as a function of the ion temperature 
gradient across the radial domain s £ [0.3, 0.8] for the helical 
a = 0 tube (left frames) and axisymmetric (right frames) 
geometry. 
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the reduced FLR suppression justifies the higher growth 
rate instabilities in the helical case. It is useful to remark 
that the g^ metric component, which does not enter 
the linear gyrokinetic equation if 9q = 0, determines the 
magnitude of g aOL via the relation g aa = (B 2 +g^ a2 )/g^ 
(which is a direct consequence of the Clebsch represen¬ 
tation of the magnetic field B = X7i[) x Va). Where the 
ITG mode peaks, z £ [— n, +7r], we have g^ a ~ 0, so that 
g aa oc l/gW. 

The different weight of the FLR suppression is not the 
only mechanism acting differently on the two geometries 
under consideration. As already mentioned in the pre¬ 
vious section, the connection length for an axisymmetric 
RFP is L c = 2iTR{q 2 + e 2 ) 1 / 2 . Such a low connection 
length is responsible for the mode suppression via Lan¬ 
dau damping at low wavenumbers, and for the high ITG 
stability threshold in axisymmetric geometry^. In a heli¬ 
cal RFP, L c is not the length of a field line corresponding 
to a poloidal turn. Contrary to the stellarator case2£, 
where L c is usually considered as the distance between 
the (stabilizing) spikes of si oc , in the RFP case the local 
shear has a more oscillatory behaviour (see Fig. [3]), so 
it makes sense to consider for L c the distance along the 
field line between two consecutive (destabilizing) peaks 
of g^. As such, we typically have L^ el L™, thus 
fcjj 1 ® 1 <C fcj[ xl . This local property makes the mode sup¬ 
pression via Landau damping less effective in the helical 
case. 

It is known that the (global) magnetic shear s influ¬ 
ences ITG mode stability. Since the q profile is essen¬ 
tially flat in the helical case for s < 0.5, this provides 
mode destabilization in the inner region. However, this 
effect is only marginal with respect to the others men¬ 
tioned above: forcing axisymmetry with a helical q profile 
yields a slightly higher growth rate with respect to the ax¬ 
isymmetric g, without modifying the wavenumber range 
where destabilization occurs, k y pi < 1, see Fig. |4|b). 

We can now compare the two helical tubes for linear 



FIG. 8. Ion heat flux (normalized to gyro-Bohm units) versus 
time for the helical geometry, on the s = 0.5 surface, for 
different ion temperature gradients, for the a — 0 and a = 
7t/ 7 tubes. 


ITG stability. Starting from the 7 (k y ) spectrum of the 
a = 0 tube, we can fix the binormal wavenumber such 
that the growth rate is maximum, k y pi ~ 1.2 in our case. 
Varying 9q (or, equivalently, k x ), it is confirmed that the 
maximum 7 corresponds to 9q = 0 (i.e., k x = 0), see 
Fig. [5}a. For the a = n/7 tube, 7 is maximum for a 
ballooning angle 9q ~ 6.27T, which corresponds to the 
maximum of g^ along the field line, Fig. [5]-a/b. In¬ 
terestingly, the highest growth rate is approximately the 
same in the two flux tubes. The dependence of k± on 2 in 
the two flux tubes is shown in Fig.[5j-c, where the value of 
6*o (or, equivalently, of k x ) is set such that 7 (^ 0 ) is maxi¬ 
mum. For the two tubes the lowest k± is approximately 
the same. For both tubes of the helical configuration a 
smaller FLR suppression is observed compared to the ax¬ 
isymmetric case. For these parameters, the structure of 
the normalized electrostatic potential as a function of the 
parallel coordinate z is shown in frame (d). 

The ITG mode growth rate as a function of the bal¬ 
looning angle 9q is shown for several neighboring flux 
tubes in Fig. [Gl It is evident that the modes have a “heli¬ 
cal ballooning” structure, being large in correspondence 
to the peaks of |Vi/>|, i.e., along the external ridge of the 
helical structure. Heuristically, the compression of the 
magnetic surfaces enhances the local temperature gra¬ 
dients, with a consequent growing instability. This is a 
common feature to any helical plasma configuration, not 
only the RFP. 

The last piece of the linear analysis on ITG stability 
refers to the radial dependence of the growth rates, Fig. [7] 
We perform this study for s £ [0.3, 0.8], this interval fully 
including the region where the transport barrier usually 
arises. Both geometries show a larger growth rate in the 
core than in the edge. While in the axisymmetric case 
the growth rate monotonically increases as s — > 0, in the 
helical case it is larger for s ~ 0.4 — 0.5, i.e., in the prox¬ 
imity of the maximum of q. The radial domain of largest 
growth rate is characterized, in both geometries, by an 
overall increased g thus by a decreased g aa and |fcj_|, 
and by a decreased global magnetic shear |s| (see Fig. 4 



FIG. 9. Parallel structure of the absolute value of the time- 
averaged turbulent potential fluctuations for the helical tube 
with o = 0 on the s = 0.5 surface. 
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of Ref. |H for an analogous study with GS2). Again, with 
respect to the axisymmetric case, the helical one clearly 
yields higher growth rates in the core, where the two 
sets of geometric coefficients largely deviate. Conversely, 
at s > 0.7 the two configurations become equivalent in 
terms of ITG stability. The footprint of the helical core 
is rapidly vanishing in the outer region: the peaks in 
g^ are no longer present and the curvature starts to 
have a tokamak-like oscillating behaviour. ITG modes 
are strongly stabilized by the (negative) magnetic shear, 
which is rapidly increasing in absolute value towards the 
edge. 


IV. ITG TURBULENCE 

Having completed the analysis of the ITG mode sta¬ 
bility for the helical and axisymmetric RFPs, we turn to 
the effect of the geometry on the turbulence. It is well 
known that, apart from the strength of the instability 
itself, ITG turbulence is largely regulated by the zonal 
flows. Therefore, again focusing on the different geomet¬ 
ric features, in this section we will also provide results 
on the linear zonal flow response for each configuration 
under study. 

We provide some details concerning the setup of the 
nonlinear simulations. For the a = 0 tube on the 
s = 0.5 surface (where q = 0.12 and s = —0.22), the 
simulations are made with the following discretization 
in the 5-dimensional space: n x x n y x n z x n v ^ x = 
128x64x512x32x8, with A k y pi = 0.05, Ak x pi = 0.055, 
and z £ [— 87 t, +87t]. For the a = 7r/7 tube, the z domain 
must be larger, so as to fully include the peaks in g xx , 
z £ [—107T,+107r], see Fig. [5] Accordingly, in order to 
allow for a full k x spectrum resolution (not centered at 
k x = 0) a larger number of x grid points is used, namely 
n x = 256, with A k x pi = 0.098. Convergence tests have 
shown that this phase space resolution is adequate. As 



1/L r . 

FIG. 10. Ion heat flux scaling versus the ion temperature 
gradient, on the s = 0.5 surface, for the helical and the ax¬ 
isymmetric tubes. The linear growth rates are also plotted 
(red), suggesting the usual nonlinear shift, which is particu¬ 
larly pronounced in the axisymmetric case. 


one can see, the most evident difference with a “stan¬ 
dard” (axisymmetric) discretization for ITG turbulence 
simulations is the large box size in the parallel direction: 
this choice is of central importance, so as to capture all 
the geometric details, especially for the a = 7 t/7 tube. 
Due to the length of the z domain and to the vanishing 
collisionality, numerical dissipations providing damping 
in the parallel direction and mimicking diffusion in veloc¬ 
ity space have been included by means of hyper-diffusion 
terms in the (z,i>||) subspace. 

In the axisymmetric case the numerical setup is less 
demanding. For the box size on the s = 0.5 surface 
(where q = 0.125 and s = —0.79) we set n x x n y x n z x 
n v n x ri/j, = 96 x 64 x 32 x 32 x 8, with A k y pi = 0.05 = 
A k x pi, and z £ [—7r,+7r]. Of course the box size, both 
in the helical and in the axisymmetric case, may differ 
from surface to surface, depending on the value of the 
magnetic shear s. 

With this simulation setup the time history of the ion 
heat flux q t for the helical geometry is shown in Fig. [8] 
for several ion temperature gradients. For two of them 
(1/= 4 and 5) we show the time trace for both the 
a = 0 and the a = n/7 tube. For the two tubes the 
same level of saturated turbulence is observed. As a con¬ 
sequence, the two tubes, as in the linear case, can be 
considered equivalent in terms of physical results, the 
ion heat flux in this case. However we remark that all 
the physical fluctuating quantities, e.g. the electrostatic 
potential cj>, the density n, the parallel and perpendicu¬ 
lar temperature Tj| and T±, exhibit a clear modulation in 
the z direction. For instance, the electrostatic potential 
is peaked where g xx peaks, see Fig. O 

In this light, from now on we can restrict our study to 
the a = 0 flux tube. Performing simulations for several 
ion temperature gradients allows us to compare linear 
and nonlinear trends, and in particular the temperature 
gradient threshold in the two cases. In Fig. [TO] we focus 
again on the half-flux surface s = 0.5. While the linear 
growth rate (dotted lines) already suggests a sharp dis¬ 
tinction between the two geometries, in a nonlinear envi- 



FIG. 11. Spectra of the time-averaged ion heat flux in k y /k x 
space, on the s = 0.5 surface, for the helical a = 0 and the 
axisymmetric tubes. 
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ronment their difference is even more pronounced (solid 
lines). In both cases a Dimits shiffeSZ exists, and the ax- 
isymmetric configuration requires ion temperature gra¬ 
dients which must be 2 — 3 times larger than the linear 
threshold for an effective ITG turbulence to take place. 
A less pronounced upshift occurs in the helical case. As 
already mentioned, these gradients in RFX-mod can be 
compared, as a reference, to the electron temperature 
profiles only, which are characterized by 1 / Lt s ~ 2 — 3 
in the core of the helical geometry, and lower in the ax- 
isymmetric one. 

Our conclusion so far is that, for an axisymmetric RFP, 
ITG turbulence is not likely to be a concern from the 
point of view of particle/heat transport, at least away 
from the edge. High current helical states, on the con¬ 
trary, are more prone to ITG instabilities and present 
a high level of ion heat flux. Besides the difference in 
the linear spectrum, discussed in the previous section, in 
Fig. [TT] we show the wavenumber spectra of the ion heat 
flux in the k y /k x space, averaged over the remaining co¬ 
ordinates and time. The helical and axisymmetric 1/Lj^ 
values are chosen to provide approximately the same ion 
heat flux, < 7 ^ ~ 1 — 2 , see Fig. [TU] In analogy with the 
linear result, the helical k y spectrum peaks at smaller 
scales than the axisymmetric one, a feature which, based 
on quasi-linear estimates, is however not enough to com¬ 
pensate for the much larger growth rates. 

We turn our attention to the effect of zonal flows 
(ZF)S£, whose role in controlling ITG turbulence in toka- 
maks is already established. In addition, in the case 
of non-axisymmetric configurations, it has been recently 



FIG. 12. Linear zonal flow response for the helical a = 0 tube 
and the axisymmetric one for k x pi = 0.2-7T, on the s = 0.5 
surface. 


showed, that although the nature of ZF can be entirely 
different compared to tokamaks, as well as among differ¬ 
ent stellarator designs, their impact is still measurable 
and beneficial. In Fig. [12] we calculate the linear ZF 
response for the RFX configuration having selected the 
normalized radial wavenumber equal to k x pi = 0.27T. In¬ 
terestingly, oscillations similar to the ones seen in the 
other stellarators are also here observed, which are at¬ 
tributed to the radial drift of locally trapped particles, 
and the residual level takes relatively small values. In 
this situation, it is suggested that the regulation of ITG 
turbulence takes place during the dynamical state of ZF 
evolution and is mildly, if at all, affected by the residual 
level. In the same figure, we show for comparison the lin¬ 
ear ZF response for the axisymmetric RFP surface. Due 
to the small safety factor and to the toroidal symmetry of 
the system, the residual turns out to be unusually large. 
Its value is in agreement with the Rosenbluth-Hinton esti¬ 
mate /Crh = 1/(1 + 1.6q 2 /e 1 / 2 ) for the undamped zonal 
flo w 30 i 31 , which for the case shown in the figure yields 
Arh|s=o.5 — 0.92. 

Coming back to the nonlinear simulations, we conclude 
with the snapshots of the electrostatic potential and of 
the density fluctuations in the x/y space in Fig. [13] for 
the two geometries under consideration (the respective 
1/LTi gradients are those of Fig. [TT]). The colour contours 
confirm the enhanced generation of the zonal flows in the 
axisymmetric configuration. 


V. CONCLUSIONS 

While a helical core has generally beneficial conse¬ 
quences on the RFP plasma performance, we have shown 



FIG. 13. Snapshot of the perturbed electrostatic potential 
and density, for the helical and axisymmetric tubes. 








































that it imposes an unfavorable effect in terms of ITG sta¬ 
bility and turbulent transport. Focusing on the role of 
geometry, thus simplifying the physical problem (we have 
assumed adiabatic electron response, vanishing plasma j3 
and collisionality, flat density profile), it turns out that 
ITG modes are localized in the proximity of the peaks 
of |V-0|, i.e. where the magnetic surface compression 
is higher. Here the local temperature gradients become 
larger, with a consequent growing instability. From the 
point of view of turbulence, the maxima of |V^| reflect 
in local maxima of turbulent fluctuations, though tur¬ 
bulence turns out to be quite uniformly spread along 
the flux tubes. With respect to the linear threshold, a 
moderate nonlinear upshift occurs, with the ion heat flux 
rapidly increasing with the ion temperature gradient and 
a consequent high stiffness of the temperature profile. On 
the contrary, for an axisymmetric RFP, already more lin¬ 
early stable to ITG modes, the nonlinear shift is much 
larger. The different behaviour is confirmed by the very 
high zonal flow residual level in the axisymmetric case, 
while helical states exhibit a lower zonal flow residual 
and a stellarator-like oscillating behaviour of the zonal 
potential. ITG turbulence is consequently expected to be 
a scarce contributor to particle transport in an axisym¬ 
metric RFP, at least without impurities. The inclusion of 
more realistic parameters of the plasma and a comparison 
with experiment are the steps to face henceforth. 
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